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Introduction 

Proper design of inlet control systems requires 
thorough knowledge of the dynamic characteristics of the 
inlet and, in particular, the dynamic response of the ter- 
minal shock to both internal and external disturbances. 
The inlet dynamic characteristics can be obtained either 
experimentally or bv means of analytical simulation. 

Experimentally, data from dynamic studies of 
mixed-compression and external compression inlets 
have been, in general, limited to disturbance frequencies 
below 50 Hz. However, in a recent investigation,^) the 
disturbance frequency was extended to 200 Hz. The re- 
sults of this investigation revealed that a resonant con- 
dition occurred in the dynamic response of this inlet at 
about 90 Hz. The dynamic characteristics were obtained 
only for a cold flow model with large internal volume. 
Ideally, the dynamic response of an inlet coupled to an 
operating turbojet engine should also be experimentally 
determined. 

A test was subsequently conducted in the Lewis 10- 
by-10 Foot Supersonic Wind Tunnel and experimental 
data were obtained for an isolated inlet coupled to short 
and long cold -flow pipe configurations and also for the 
inlet coupled to an operating turbojet engine. The iso- 
lated inlet and the inlet-engine configurations were sub- 
jected to an internal disturbance at frequencies up to 
240 Hz and also to an external disturbance at frequencies 
up to 15 Hz. Data were also obtained for the inlet with 
and without boundary layer bleed in the throat region. 

A mathematical description of the inlet dynamic be- 
havior is also briefly presented. The analysis combines 
a linearized set of equations across th shock with a 
relatively exact solution of the subsonic duct wave equa- 
tions. The solutions for the wave equations are written 
in terms of time shifts and then expressed in the fre- 
quency domain. Data from a digital computer solution of 
the equations for the inlet frequency response are pre- 
sented. The complete analysis is presented in Ref. 2. 

A comparison between the experimental and analytical 
results is made. 

Apparatus and Procedure 

The experimental investigation was conducted in the 
Lewis 10-by-10 Foot Supersonic Wind Tunnel at a free 
stream Mach number of 2.5. An isometric view of the 
inlet is shown in Fig. 1(a). The inlet was designed for 
Mach number 2. 5 operation with an inlet cowl lip diam- 
eter of 0. 465 meter and with 60 percent of the total 
supersonic area contraction ratio occurring internally. 
The inlet matched the airflow requirements of a J-85-13 


turbojet engine. Boundary layer bleed on the cowl and 
centerbody were provided in the inlet throat region. A 
detailed description of the inlet steady state aerodynamic 
performance is reported in Ref. 3. 

Internal and external disturbance devices were used 
to perturb the inlet flow. The internal disturbance de- 
vices generated symmetrical and unsymmetrical flow 
perturbations. The symmetrical disturbance was ob- 
tained by sinusoidally oscillating three of the six inlet 
bypass doors at frequencies up to 140 Hz. The locations 
of the hydraulically operated bypass doors are shown on 
Figs. 1(a) and (b). The unsymmetrical disturbance was 
obtained from a rotary valve positioned in place of one 
bypass door, Fig. 1(b). The valve exit area was de- 
signed to vary sinusoidally during constant speed rota- 
tion. The external disturbance was a trapezoidal gust 
generator plate positioned above and upstream of the in- 
let and oscillated through 1° of arc at frequencies up to 
15 Hz, Fig. 1(b). Construction and operation of the ro- 
tary valve and gust generator plate is described in Ref. 1 , 

The inlet shock position dynamics were obtained 
from eight static pressure transducers located in the in- 
let throat, Fig. 1(a). The analog signal of each trans- 
ducer indicated the time of shock crossing. The times 
at which the shock crossed the various transducer posi- 
tions was then used with the method of least mean square 
to obtain amplitude and phase shift of the shock position 
dynamics^ 1 ). The response of the throat exit static 
pressure was recorded on analog tape by a transducer 
shown in Fig. 1(a). Calibration of the tube -transducer 
assembly allowed for correction of the pressure fluctu- 
ation signals due to tube assembly dynamics. 

Experimental data for the external and internal dis- 
turbances were obtained for the inlet coupled to short 
and long cold pipe configurations and also for the inlet 
coupled to an operating turbojet engine. The short and 
long cold pipe configurations had the internal flow choked 
at the diffuser exit and 2.36 meters further downstream 
at the model exit station, respectively. In addition to 
these configurations, dynamic data were also obtained 
with and without boundary layer bleed in the inlet throat 
region . 

The dynamic responses presented in the figures are 
in terms of amplitude ratio and phase shift. The ampli- 
tude ratio is defined as the amplitude at the test fre- 
quency divided by the amplitude at a frequency of 1 Hz. 

Faired curves rather than data points are presented 
on Figs. 2, 4, 5, and 6 so that the comparisons between 
the different inlet configurations can be made easily. 

The data points are presented on Figs. 3, 9, 10, and 11. 
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Experimental Results 
Internal Disturbances 

The dynamic responses of shock position to a sym- 
metrical internal disturbance for the three inlet config- 
urations are presented in Fig. 2. For these configura- 
tions, the inlet throat region had normal boundary layer 
bleed on the cowl and centerbody. The ratio of the 
change in shock position to the change in weight flow 
(AX g /Aw) was about 0.122 meter per kilogram per 
second at the reference frequency of 1 Hz. The distur- 
bance amplitude was ±0.40 kg/sec for the cold flow con- 
figurations and ±0.30 kg/sec for the inlet-engine config- 
uration. The maximum capture inlet flow rate was 
12.90 kg/sec. 

All three configurations exhibited a resonant condi- 
tion at about 50 to 55 Hz. Little or no change in reso- 
nant frequency was noted when the location of the inter- 
nal flow choke point was changed and when the engine was 
coupled to the inlet. The shock position dynamics for 
the choked diffuser exit configuration and the inlet- 
engine configuration showed minor differences except in 
the vicinity of 30 Hz. Increasing the diffuser volume 
decreased the magnitude of the amplitude ratio more 
rapidly at lower frequencies and accentuated the reso- 
nant condition as indicated by the solid curve. First- 
order lumped parameter analyses, such as reported in 
Ref. 4, predict the data to about 25 or 30 Hz. However, 
this type of analysis cannot predict resonances such as 
those observed in this test. 

The variation in phase shift with disturbance fre- 
quency appeared to be similar for the choked diffuser 
exit and the inlet-engine configuration. The phase shift 
for the configuration with the model exit choked de- 
creased more rapidly at low frequencies but a leading 
effect is noticed just prior to the resonant frequency. At 
the higher frequencies the phase shifts appear to be simi- 
lar for all three configurations. 

The dynamic response of shock position to an unsym- 
metrical internal disturbance is presented for fre- 
quencies up to 240 Hz, Fig. 3. The internal flow was 
choked at the model exit for this configuration. The data 
presented here are for amplitude ratio only and were ob- 
tained by the inlet unstart method described in Ref. 1 . 

The incremental change in weight flow through the rotary 
valve was constant at about ±0.13 kg/sec over the entire 
frequency range. 

The location of the shock position at unstart is up- 
stream of its normal location in the inlet throat where 
the usual shock position response data are taken. This 
difference plus the unsymmetrical nature of the internal 
disturbance could cause the differences in response noted 
in comparison of this data to that presented in Fig. 2 for 
the same configuration. Resonant conditions appear at 
about 60, 125, and 210 Hz. The magnitude of the ampli- 
tude ratio at the resonant condition of 60 Hz was equal to 
the magnitude at 1 Hz. 

Although the resonant conditions were at frequencies 
greater than 40 Hz, the test model under consideration 
was of moderate size (0. 465 meter cowl lip diameter) 


and therefore the resonant frequencies should be reduced 
by a scale factor for comparison with larger models. By 
limiting the permissible gain for a stable control system, 
the presence of these resonances can affect controlled 
system performance at relatively low frequencies. 

Figure 4 presents a comparison of the dynamic re- 
sponses of shock position and throat exit static pressure 
for a downstream symmetrical internal disturbance in 
the inlet-engine configuration. At the reference fre- 
quency of 1 Hz, the change in static pressure at the 
throat exit was 6.89x10 s newton per square meter for the 
corresponding change of 0.59 kg/sec in the bypass weight 
flow. 

In general, the dynamic response of static pressure 
just downstream of the inlet throat should lead the shock 
position dynamics in a first order manner. Similar val- 
ues in amplitude ratio were noted for the two response 
curves to a frequency of about 70 Hz, and the two respon- 
ses are within a factor of 0.707 to frequencies beyond 
90 Hz indicating a first order corner frequency larger 
than 90 Hz. 

Comparison of the two response curves in phase 
shift (dashed and solid curves) shows a larger difference 
than would be expected from a first-order relationship 
whose corner frequency is larger than 90 Hz. However, 
consideration of the dead time between the average shock 
location in the throat and the throat exit static tap (15 cm 
downstream) shows that the difference in phase shift of 
the two responses (dashed and broken curves) to be less 
than 45° for frequencies beyond 90 Hz. The analysis 
presented later in this report predicts the corner fre- 
quency to be about 120 Hz. 

The effect of inlet throat bleed on shock position dy- 
namics for the inlet-engine configuration is presented in 
Fig. 5. The resonant frequency was reduced from 55 Hz 
to approximately 40 Hz when the downstream throat bleed 
was sealed. The ratio of the change in shock position to 
change in bypass weight flow (AX s /Aw) for the inlet with 
the downstream throat bleed sealed was 4. 4 times the 
(AX s /Aw) for the inlet with full throat bleed at 1 Hz. 
Therefore, the terminal shock of the inlet with sealed 
bleed was more sensitive to small downstream distur- 
bances while for the inlet with the full throat bleed, the 
terminal shock was less sensitive but subject to a more 
pronounced resonant condition. Small differences in the 
phase shifts of the two configurations were noted with the 
phase of the sealed bleed configuration laging that of the 
full bleed configuration. 

External Disturbance 

The response of shock position to an external distur- 
bance for the inlet-engine configuration with the down- 
stream throat bleed sealed is presented in Fig. 6. The 
ratio of the change in shock position to change in gust 
plate angle of attack (AX s /A<p) was 0.070 meter per 
degree. The external disturbance perturbed the flow in 
both the supersonic and subsonic portions of the duct. 

The gust plate produced simultaneous changes in free 
stream Mach number and flow angle of attack in the flow 
field ahead of the inlet. 
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The data of Fig. 6 indicates a moderate reduction 
in the amplitude ratio to about 10 Hz. From 10 to 15 Hz 
the response in amplitude was flat or slightly increasing 
in magnitude. This would suggest that a possible lead 
component in the amplitude ratio may have been present 
after 15 Hz. Little difference in phase shift was noted 
for the two test conditions over the range of test fre- 
quencies. In Ref. 1, inlet response data resulting from 
an external disturbance indicated a lead characteristic 
in both amplitude ratio and phase shift causing the inlet 
to be sensitive to higher frequency external disturb- 
ances. 


Analytical Results 

The one dimensional mathematical analysis of shock 
position dynamics in a mixed compression inlet can be 
conveniently subdivided into two problems. First the 
interaction of shock position with adjacent variables 
and, secondly, the subsonic duct analysis. Figure 7 
shows the idealized inlet considered in the analysis. 

The inlet was treated mathematically as a duct with 
varying area followed by a duct with constant area. The 
upstream end of the duct is bounded by a moving normal 
shock and the downstream end by a fixed choke station. 

The analytical results here presented are based on 
a linearized set of equations across the normal shock 
combined with linearized wave equations for the sub- 
sonic duct. The linearized wave equations can be inte- 
grated only over a constant area section. The diverging 
portion of the subsonic diffuser is thus mathematically 
represented as three constant area sections. 

Previous dynamic analyses of inlets have applied 
various techniques to the problem. Reference 5 pre- 
sents a detailed analysis of the lumped volume Helm- 
holtz resonator concept. This representation, however, 
does not predict the repeated resonances typical of a 
distributed parameter system. The analysis applied in 
this paper does predict such resonances without being 
significantly more complicated in application than the 
Helmholtz -resonator concept. Efforts have also been 
directed toward the application of the method of char- 
acteristics to the problem. This technique should pro- 
vide accurate results but with a significant increase in 
computation time and analytical complexity. 

The linearized equations used in this analysis can 
provide a closed form solution for each section and 
when combined yield an easily calculated frequency re- 
sponse. It is also possible to program the equations in 
a simulation language such as Mimic or on an analog 
computer for direct control system simulation. 

The equations for both shock motion and the sub- 
sonic duct are written with flow, total pressure, and 
entropy as the dependent system variables. These pa- 
rameters have the advantage of being constant in the 
steady state, and thus considerably simplify the small 
perturbation analysis. A more thorough development of 
the analysis is presented in Ref. 2. Mass, momentum, 
and energy balances taken across a normal shock can be 
modified to account for shock motion. With the assump- 
tion of perfect gas relations and negligible changes in 


duct area acorss the shock, a linearized set of equations 
is presented. Restricting consideration to disturbances 
originating downstream of the normal shock, the dynam- 
ic relations of equations (1) are derived: 

?2 = ~ C l[ L T a' + C 2 S ] X S 1 


w 2 = C 3 sX s 

s 2 = Ci [l t £ -c 3 s]5c s 


( 1 ) 


where 


and 


Ci = 


(7M2 - 1) 


5<1 +0.2M 2 ) 1 / 2 T 
C 2= ' Lt 


C 3 " a T Mj 


(Mi - 1) 


(1 + 0. 2Mp 


1/2 T 


and a ~ over a quantity designates a small perturbation 
variable, that is, P ~ AP T /P T , w = Aw/w, and S = 
AS/R. 


In Ref. 2 the conventional compressible flow conti- 
nuity, momentum, and energy equations for quasi-one 
dimensional flow are formulated as wave equations in 
terms of the chosen system variables total pressure, 
flow, and entropy. With the assumption of invicid, non- 
heat conducting flow of an ideal gas it is shown that: 


g 5+p T + _1_ + aM(M + 2) 6+S _ Q 

apj 6t A p 6t 7R 5t 
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The above equations are then linearized for small per- 
turbations about a steady flow condition and with the sub- 
stitution of: 


a - ■ 7M -. 

(M 2 + 5) 

p M(M + 2) 
<M 2 + 5) 

and 

M(2 - M) 

7 

<M Z + 5) 
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( 8 ) 


equations (2) become: 


A + ( V + a) — 

fit fix 


-A + (v - a) -A 

fit fix 


[jP T + aW + /3sJ = 0 

[ 


P T + aw - ysj = 0 


[5"Er-' 


(3) 


The solutions of equations (3) may be written in terms of 
the Laplace transform operator s as: 


r~ -v ~ i 

= e- (TS 

L ~ 

P . - + Qfw + /3S - 

L n+1 n+1 n+l_ 

_ P n + mv n + /3S n| 

[ P n+1 - ""'n+1 - '>' S n+l 

^e 4 * 78 

p n - - ys n ] 


y = e_0s 5i 


where a = L/(v + a), r = L/(a - v), and 0 = L/v and 
L is the linear distance between stations n and (n + 1). 
The integration performed between equations (3) and (4) 
is only valid for a constant area duct. The variable area 
portion of the subsonic diffuser is assumed to be made 
up of three constant area sections. The boundary con- 
ditions between the sections pose no problem since total 
pressure, flow and entropy may be regarded as constant 
across the interface. 


Equations (1) and (4) when combined with appropri- 
ate terminal conditions for the subsonic duct and distur- 
bance equations written in terms of the selected vari- 
ables, provide a mathematical model which can be ap- 
plied to a mixed compression inlet. 

In the mathematical analysis of Ref. 2, the effects 
of performance bleed were not considered. For the test 
data here presented, however, it is necessary to include 
the performance bleed. If perturbations in upstream 
flow are included in the second of equations (1), it be- 
comes: 


Aw Aw. AX 

w w T 


(5) 


Figure 8 shows a normal shock standing in a region with 
performance bleeds. In this figure stations 1 and 2 
are just beyond the range of normal shock travel, while 
stations la and 2a are adjacent to the normal shock 
in the bleed region. With the nomenclature of this fig- 
ure, equation (5) becomes: 


Aw 2 _ Aw 2a A ^2b 

W w w 

If the effects of small perturbations in pressure are ne- 
glected, the bleed flows are proportional to shock posi- 
tion. 

- AX S 

w ib = w ib + %b TT 
T 


AX, 

Aw ib = K ib T 


s 


and 


Aw 2b = -K. 


T 

AX, 


2b ' 


s 


(9) 

( 10 ) 


Substituting into equation (6) results in 


Aw 2 

w 


Awi 

w 


= < K 2b- K lb>-[^ +C 3 s 



or neglecting upstream disturbances 

w = <C B + C 3 s)X s (11) 

Steady state test data obtained for the performance 
bleed flow as a function of shock position indicates the 
linear relation predicted by equation (11). 


In the actual application of the analytical technique 
to a real problem, the principal difficulty arises in 
determining the appropriate value of A'/A (i.e. , 1/A 
dA/dx in the vicinity of the normal shock). Because of 
boundary layer effects, the actual geometric value does 
not give satisfactory results. It was found empirically 
that A 1 /A = 0. 013/cm gave good agreement for the low 
frequency gain and fair frequency response agreement. 
This value was thus used in all of the analytical data 
presented. The bleed coefficient (0.001 51 /cm) was 
based on actual experimental data. Although the curves 
presented show fair correlation between experiment and 
theory, the results could be improved by better selec- 
tion of these two parameters. 


Figure 9 presents shock position frequency response 
obtained for the experimental inlet with the long cold 
pipe attached. Both the test data and the analytical re- 
sults exhibit the repeated resonances typical of a dis- 
tributed parameter system. 


Aw 2a . Aw la 
w w 


AXg 

+ c 3 s ~d 


( 6 ) 


The flow at station 1 is related to that at station la by 


or 


w = w, +w lL 
1 la lb 


AWi -AWjq | Aw lb 
WWW 


(7) 


where w-^ is the bleed flow between stations 1 and 
la. Similarly: 


Figures 10 and 11 present frequency response data 
for shock position and static pressure just downstream 
of the shock. These two figures are also for the exper- 
imental inlet, but with the cold pipe replaced by a 
choked station at the end of the variable area portion of 
the duct. In both plots the phase agreement is good but 
although the analytical results exhibit a plateau in the 
amplitude ratio plots they fail to exhibit the distinct 
resonance of the test data. One possible source of the 
discrepancy arises from the fact that the actual disturb- 
ance device is separated from the inlet duct by a vol- 
ume. This volume has not been included in the analyti- 
cal model. 
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Concluding Remarks 


shock position 


Dynamic response of the shock position resulting 
from a symmetrical internal disturbance showed that: 


1. Resonance was observed at about 55 Hz. 

2. Little or no change in the resonant frequency was 
observed when the internal flow choke points were 
changed or when the engine was coupled to the inlet. 

3. The dynamic response of the shock position was 
observed to be similar for the compressor face choke 
point and inlet-engine configurations. 

4. The additional volume of the inlet configuration 
with the internal flow choked at the model exit resulted 
in a more rapid decrease in amplitude ratio. 

5. Sealing the downstream throat bleed resulted in 
increased sensitivity of the terminal shock to down- 
stream disturbances. 

6. Dynamic response of the shock position to an ex- 
ternal disturbance was nonresonant over the test fre- 
quency range for the inlet-engine configuration. 

7 . Comparison of the analytical and the experimen- 
tal data Indicates that fair agreement can be achieved 
between the two. Further work should include the ef- 
fects of the volume associated with the disturbance de- 
vice. In addition, a finer lumping of the diverging por- 
tion of the subsonic diffuser should be investigated to 
determine the effects on accuracy. 


Symbols 


A area 


A* dA/dx 


a 

C B 

Ci,C 2 ,C 3 

g 

K lb> K 2b 

L 

M 

n 

P 

R 

S 

s 

t 

v 

w 


speed of sound 

bleed flow coefficient 

coefficients with appropriate dimensions 

acceleration due to gravity 

bleed flow coefficients 

length 

Mach number 
station number 
pressure 

universal gas constant 
entropy 

Laplace transform operator 

time 

velocity 

flow rate 


X s 

x space coordinate 

ce,/3,y functions of Mach number defined in text 

p density 

cr, r, 6 delaytimes 

ip gust plate angle of attack 

The subscript T denotes total conditions. A bar over 
a variable denotes the steady state value of the variable. 
The symbol A preceding a variable denotes a small 
perturbation quantity. In general, numerical subscripts 
denote the station associated with the variable as defined 
in Fig. 8. 
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Figure 2. - Dynamic response of shock position to sym- 
metrical internal disturbance. 
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Figure 3. - Dynamic response of shock position to unsymmet- 
rical internal disturbance. Model exit choked; inlet unstart 
method. 
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Figure 4. - Response of shock position and throat exit static 
pressure; inlet-engine configuration. 



Figure 5. - Effect of inlet throat bleed on shock position 
dynamics; inlet-engine configuration. 
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Figure 6. - Shock position response to external disturbance. 
Inlet-engine configuration; throat bleed sealed. 
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